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Abstract 

We propose a method for the approximate computation of the Green function of a 
scalar massless field f subjected to potential barriers of given size and shape in spacetime. 
This technique is applied to the case of a 3D gaussian ellipsoid-like barrier, placed on the 
axis between two pointlike sources of the field. Instead of the Green function we compute 
its temporal integral, that gives the static potential energy of the interaction of the two 
sources. Such interaction takes place in part by tunneling of the quanta of cj) across the 
barrier. We evaluate numerically the correction to the potential in dependence on the 
barrier size and on the barrier-sources distance. 
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I. Introduction. 


In Quantum Field Theory it is useful in several occasions to have a general expression for the 
Euclidean two-point correlation function of a massless scalar held (j) in the presence of potential 
“barriers” in spacetime of the form 





( 1 ) 


Jn{x) being the characteristic function of the 4-region hi where the potential has support (= 1 
for xeQ, Jq = 0 elsewhere). The region hi can be multiple connected, thus representing several 
barriers placed at different points in spacetime. 

Possible applications are connected for instance to the fact that a potential of the form 
(|l]) represents a localized imaginary mass term (m^ < 0) in the action of the scalar held (j). 
Terms of this kind can be present in cosmological models with inhationary helds. It is also 
known that every quantum held with non-vanishing vacuum expectation value (VEV) has a 
global imaginary mass term in its lagrangian |l|], which couples to the gravitational held as 
a cosmological term; one can show ^ that if the VEV is not constant but depends on x, it 
becomes a local cosmological term for the gravitational held. 

More generally, suppose we have a system of two interacting helds and regard one of them 
(or its VEV) as a hxed external source. The coupling term of the two helds becomes a local 
constraint for the dynamical held, a sort of external potential localized in the regions where the 
external held has support. It is therefore important to study the tunneling of the dynamical held 
through these regions, that is, its Green functions. Note that in systems like this translational 
invariance is generally lost. 

It is easy to check that the potential V (0) in eq. (|^) implements in fact a constraint in the 
functional integral of the held: writing this integral as 


z = 


J (00)^ — y d‘^xV(0) , 


one sees that for large ^ the square of the held is forced to take the value 0 q within the region 

Q. 

For a characteristic function Jn{x) like the one specihed above we say that the constraint is 
imposed in a “sharp” way: in spacetime the potential barrier looks like a step at the boundary 
of G. Smoothing Jq we can obtain a smooth potential barrier. In the following we shall be 
more interested in this second case. 
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Note that the potential (|^ has the shape of a double well, as long as considered only a 
function of the held 0, but regarded as a function of x it is positive and reminds much more a 
barrier. 

Let us focus on the case of weak helds, such that 0*^ can be disregarded with respect to 0^. 
If the product 7 = .^0q is small, then the effect of the barriers on held correlations is small, too, 
and can be treated as a perturbation. One can solve the equation for the modihed propagator 
G'{xi,X 2 ) = (0(a;i)0(a;2))v in closed form (see the Appendix), hnding that G' is given by a 
double inverse Fourier transform, with the direct transform of Jq evaluated at {p + k): 

G{xi,X2) = G^{xi,X 2) + jG'ixi, X 2 ); 

G\xi,X 2 ) = J 

G'{x,,X 2) = Jd^pJ ( 2 ) 

In hnite-dimensional quantum mechanics computing G\xi,X 2 ) corresponds to compute the 
Feynman transition amplitude, related in turn to the system’s wavefunction in the presence 
of barriers. In held theory the intuitive meaning of G\xi,X 2 ) is less immediate. However, we 
can derive from G'{xi,X 2 ) a quantity with a direct physical interpretation: the static potential 
t/(xi,X2) of the interaction of two pointlike sources qi and q 2 of the held 0 at rest. This 
interaction is mediated by the exchange of quanta of 0. If the barriers are placed somewhere 
between the sources, the interaction is clearly ahected, but it still takes place - provided the 
product 7 is small - with the quanta of 0 “tunneling” through the barriers (or passing over the 
wells, depending on the interpretation). 

The leading contribution to the static potential f/(xi,X2) is obtained from (|^) as follows 
[§]. First one dehnes Jn{x) as the product of a 3D function jo(x) and a function constant in 
time, then one integrates over ti and ^2, multiplies by qiq 2 and divides by —T, taking the limit 
for T —> 00. The result is 

17(Xi,X 2) = t/°(Xi,X2) +7l7'(Xi,X2) = 

= I I - 7(27r)®gig2 [ dp [ ' ^3) 

|xi—X2I J J 

This formula is easily generalized to the case of N charges qi,...,qN, placed respectively at 

Xi, ...,XAr. 

A limit case of the physical situation we are considering is represented by the electrostatic 
potential of pointlike charges in the presence of perfect conductors. In this case the held is 
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exactly zero within the region f2, and has sharp bonndaries — thus jo(x) is a step function 
and jo(p) a strongly oscillating function. Eq. (|^) could be applied to this case only if the 
parameters 0o and ^ could be chosen in such a way that 0o ^ 0 and ^ oo, the product 
7 = still being finite and small. We know, however, that usually in an electrostatic system 
the change in potential energy due to the presence of perfect conductors is not just a small 
correction. (It can be computed exactly, in principle, solving a classical field equation with 
suitable boundary conditions.) 

The case of interest here is actually more subtle. In the following jn(x) is supposed to be a 
smooth function and both 0 o and ^ are taken to be finite. The field square has only a certain 
probability to be equal to 0 q within This probability is maximum at the center of and 
decreases towards the boundary of ff. Since is smooth (a gaussian function), its Fourier 

transform jn{p) is smooth too, and the integral can be computed numerically. 

It is interesting to study U' in dependence on the geometrical features of the barrier and 
on the position of qi and q2 with respect to it. Take, for instance, a finite size barrier (gaussian 
ellipsoid, see Section II) lying on the axis joining xi to X2. We may expect that if one of the 
two charges is close to ff, then |f/'/[/°| is larger, decreasing if both charges are far away from 
— or if is not on their axis. This behavior is confirmed and specified by our numerical 
results. 

The paper is organized as follows. In Section II we compute the leading correction to the 
static potential for a barrier with the shape of an ellipsoid. Due to the peculiar behavior of the 
integrand, the procedure for numerical integration is not trivial and requires some care. We 
describe it in detail. Results are given in Section III. They concern in particular the dependence 
of the correction to f/(xi,X2) on the geometrical setting (size of the barrier and its position 
with respect to the pointlike sources). Far from exploring all the conceivable variations and 
related phenomenology, the main aim of this work is to show that the general technique can be 
successfully applied to real cases. 


II. The case of two static sources. 

Let us focus now on a configuration with two static sources and one barrier only. We choose 
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our reference frame in such a way that the sources he on the z-axis: 


Xi = ( 0 , 0 ,Li); X 2 = ( 0 , 0 ,-La). 


The spatial shape and size of the barrier are dehned by the function 


jn(x) = exp 




( 4 ) 


This means that the region is like an ellipsoid centered at the origin, with simmetry axis 
along Oz, radius of the order of a and thickness of the order of h. We suppose that a > b, thus 
the ellipsoid is “squeezed” on the x|/-plane. More precisely, the region itself is not sharply 
dehned, but the surfaces where jn(x) is constant are ellipsoids. For instance, on the surface 
dehned by 

x‘^ + y‘^ z‘^ 

a 2 ^ 

the function jo(x) is constant and equal to e~^. 

The Fourier transform of is 

Jo(p) = 7r3/2a26exp|^ 


« iPx + Pv)-b p: 


2^2 


( 5 ) 


The charges qi and q2 can be taken to be unitary and the distances Li and L2 expressed as 
multiples of the ellipsoid radius a: Li = riia, L2 = n2a. From ( 1 ), (|]) we obtain 

^ikznia—ipzn2a 


17'(xi,X 2) = —(27r)®7r^/^a^6 J dp Jdk- 


k 2 p 2 


■ exp 


~a^ip + k)ly - + Kf 


where ^xy denotes the component of a vector V in the plane xy. In the following we shall be 
most interested in the case with one charge far from the barrier (ni S> 1 ), while the other charge 
is close to it (typically in our numerical calculations n2 ranges between 1 and 15 ). Accordingly 
we set = e, n2 = n. After rescaling —>■ ek^ we obtain 

„ ^ikza—ipzua 


£/'(x„x,) = dpj dk 


exp 


~a^ip + k)ly - + ek,f 


Then we eliminate any further dimensional parameters by rescaling k ^ 2 k/a and p —>■ 2 p/ a, 
obtaining 

^ikz—ipzTi 


t/'(x.,x,) = 


exp 


-iP + ^)xy- P {Pz+eKY 


where p = b/a is the ratio between the thickness b and the radius a of the ellipsoid. 
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Next we introduce the polar variables 6k, dp, (pk, and (pp. In the following k and p will not 
denote four-vectors anymore, but |k| and |p|, respectively. The square of the component of the 
vector (p -|- k) in the xy plane is 

(p + k)^j^ = sin^ Op -|- sin^ 9k + 2pk sin dp sin dk cos{(pk — (pp)- 

The other components are 


kz = k cos dk, Pz = P cos dp] 


K + ky = ^ 

Finally, introducing the variables 

s = cosdk, t = cos 9p, (p = {(pk — (pp) 

one obtains, remembering that the integrand is even in s, t, the following basic formula: 

N (27r)^° ab ^ , f°° ,, , cos{ks - npt) 

U (xi,X2) =- — — 271 / d(p ds dt dk dp- - — -^ x 

^/7^ Li Jo J-i J-i Jo Jo l — s^{l — e‘‘) 


X exp 


-p‘^{pt + eksY — p^{l — t^) — — s^) — 2pk cos (p\J (1 — (1 — s^) 


(27r)^^ ab f‘^'^ 

- d<p ds dt dk dp f{(p, s,t,k,p]e, p,n) 
1 Jo J-l J-l Jo Jo 

F{e,p, n). 


yT L 
{27iY^ ab 
71 Li 


( 6 ) 


A. Preliminary study of the integrand. 

It is important to discuss in advance the case in which p and £ take values much smaller 
than 1, that is, hi is very thin and the distance of the hrst charge from hi is much larger than 
a. When t and s approach -|-1 or -1, for small values of p the integral over k and p converges 
very slowly at inhnity and the factor cos{ks — npt) performs a large number of oscillations. For 
very small £ there are many more oscillations in k than in p. (In the limit p ^ 0 the integral 
makes sense only as a distribution. We shall never approach this limit, however.) 

Let us set, for instance, s = 1, f = 1 and (p = t: 12 in the argument of the exponential in (^). 
We obtain the exponential factors 


exp 


—p^{pt + eks)"^ 

= exp 

-{pp? 

exp 

-{pek)'^ 

exp 

—2p^ekp 


(7) 
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The first factor on the r.h.s. of (J^) has a range in p of the order of p~^ and the second factor has 
a range in k of the order of {pe)~^. The third factor has a range in p, for fixed A;, of the order of 
{p^/ek)~^ and a range in k, for hxed p, of the order of {p\/sp)~^. Fortunately this latter factor 
is not relevant: if its range is larger than the other two ranges then it does not play any role; 
if it is smaller then it is sufficient to refer to the other ranges. 

As soon as and go away from 1, the number of oscillations of the integrand decreases. 
For instance, setting s = t = 0.98 we obtain the exponential factors 


exp 


— {ppY — — 2p^ekp 


exp 


-0.04-0.04 


When p is much smaller than 1 the range of this product is determined by the second exponential 
and does not depend on p. 

It is also easy to take into account the term proportional to cos 0. After setting 0 = tt that 
term gives a positive contribution to the argument of the exponential; thus studying the range 
of the resulting expression we obtain an upper limit valid for any 0. 


B. Integration domains. 

Independently of the considerations above, it is possible to plot the integrand /(0, s, A, k, p; e, p, n) 
for several different values of p and £ and check the ranges of the exponentials. In order to 
better control the oscillations of /, we study it in 4 different domains of the variables s, t: 

Domain 1: t,se [0,1 — a]; 

Domain 2: Ae [0,1 —a]; se[l —a, 1]; 

Domain 3: A e [1 — a, 1]; s e [0,1 — a]; 

Domain 4: A, s e [1 — a, 1 — a]. 

A typical value of a employed in the program is a = 0.02. The total integration domain in 
s, A is obtained by “reflecting” each of the domains above with respect to one axis and then 
reflecting again the result with respect to the origin (s —»• —s, A —>• —A, s,t —>■ —s, —A). In each 
domain i there is a maximum value for the variables k and p, beyond which / is equal to zero 
for any practical purpose. Denoting by Ki and Pi these ranges. Table 1 shows the results found 
for some considered values of £ and p. 

Since the integration over k and p is extended to wide ranges, the most reasonable technique 
for the numerical computation of the integral (^) appears to be a Montecarlo sampling of the 
integrand. The sampling algorithm evaluates the average value of / in each domain, extending 
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the values of k and p up to the maximum range necessary for that domain. At the end the 
global average is computed, weighing each single average with the ratio between the domain 
volume and the total volume. Denoting by fi the average of / in the domain i and by Vi the 
domain volume we have 
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^ = E = 8^ (1 - + «(1 - a)K2P2f2 + «(1 - a)KsP3h + a^K.PJ, . (8) 


i=l 


III. Results of the numerical integration. 

The contributions of the Domains 2, 3 and 4 to the integral F (compare (^) are found to 
be small with respect to the contribution of Domain 1. The fluctuations of the average of / in 
Domains 2 and 4 (where approaches 1) may be very large. In order to achieve a sufficient 
precision these regions have been sampled with a large number of points (up to ~ 10^°). The 
standard routine “ran2” was used for random numbers generation. 

The dependence of the integral F on the parameters e and p is very weak, thus U' depends 
on a, h and Li mainly as ah/Li (see eq. (|[)). The study of the dependence of U' on n is more 
difficult, because this dependence is entirely contained in the integral F and can be evaluated 
only numerically. One needs to insert in the program a cycle which samples the integrand for 
different values of n, typically between 1 and 15. This is possible because the ranges P*, Ki do 
not depend on n. 

The numerical evaluation of P as a function of n in the range n = 1...15, with e: = 0.1 and 
p = 0.3, gives the results shown in Fig. 1. With p = 0.1 and p = 0.032 one obtains very similar 
results, thus confirming the weak dependence on p (Fig. 2). As expected varying £ does not 
affect much the value of F either, since the dependence on the distance Li is already factorized 
out of the integral (compare Fig. 3). 

Fig.s 1, 2, 3 reveal an exponential behavior of P(n) of the form 

F{n) ~ exp(—mn + q) + h. 

It is also clear just from the graphs that the exponential decrease of P for large n leaves an 
asymptotic value F = b, with b in the interval 0.1-0.2. This is an interesting behavior, as it 
means that the “shadow” produced by the barrier in the static field of the two sources has a 
long, constant tail. 




A least-squares fit of the data gives the results of Table 2. Excluding from the fit the first 
two points {n = 1 , 2) we obtain better estimates for the distribution tail and for b. The errors 
on the parameters of the fit, in particular those on b, are small. They can be estimated knowing 
that the least-squares sum of the percentual errors S = ~ [exp(—mn + q) + 6]/F(n)}^ 

has a minimum value Smin ~ 0.05 and that its second partial derivatives at the minimum are 
of the order of 0 ~ 4 • 10 ^, 0 ~ 10 ^, 0 ~ 10 . 


IV. Conclusions. 


Our technique for the computation of the Green function and the static potential of two 
pointlike sources appears to work well for weak helds, yielding reasonable results. The method 
is based upon a double 3D Fourier transform of the function which represents size and position 
in space of the potential well or barrier. This double transform is necessary, due to the lack of 
translational invariance of the system. Its numerical evaluation requires a preliminary analytical 
study and a subdivision of the integration volume in a few domains, because the range of the 
real exponential factors in the integrand varies considerably. 

We studied the case of a smooth barrier with the form of a gaussian ellipsoid in coordinate 
and momentum space. For values of p and e not much smaller than 1 a good precision was 
obtained, (p is the ratio between the lengths a and b of the ellipsoids axes and £ is the ratio 
between the length of the major axis a and the distance Li of the hrst source from the ellipsoid.) 

Denoting by n the distance of the second source in units of the major axis, we found that 
the correction to the interaction potential along the line joining the two sources and the barrier 
has the following form (compare eq.s (|), (D): 


U = G° + 7G' = 


Xi - X2 


|Xl - X2I 


( 27 r)^^ 706 
TT Li 


F{e,p,n) 


1 - 


(27r) 


TT 


-706(1 -f ne)F{e, p, n) 


The function F depends very weakly on p and e. Its dependence on n is displayed in Fig.s 1-3 
and shows an exponential decay followed by a constant tail. 

The behavior summarized above is interesting in itself, being the result of a sort of “tun¬ 
neling” of the scalar held through a region where it is constrained or has imaginary mass. We 
have seen that the local imaginary mass term affects the propagation of the held also outside 
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the region ^2 where it has support. This feature is easily understood from the physical point of 
view; we gave here a method for its quantitative evaluation. 
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Appendix: Proof of the expressions for G', U'. 


We give here the proof of eq.s (H) and @ of the main text. Expanding the square in (|I]) we 
obtain for W[J] 


W[J] = j d[4>\ exp 


(Tx 


{d(t)f 





The last term in the square bracket is constant with respect to 0(x) and its exponential can be 
factorized out of the functional integral. In a first instance - for weak helds - we can disregard 
the term. We are then led to consider a quadratic functional integral, and the “modihed 

propagator” G{x,y) = (0(a;)0(j/)) j, which by dehnition satishes the equation 


dl + lMx) 


G{x,y) 


-( 27 r)V(x - y), 


(9) 


where 7 = > 0. Let us focus on the case when (j^l = 0 inside the regions hlj and let us take 

the limit 0o 0 and ^ cx) in such a way that 7 is hnite and very small, so that the term 
'yJn{x) in eq. (j^) constitutes only a small perturbation, compared to the kinetic term. Then 
we can set 

G{x, y) = G^{x, y) + 7G'(x, y), 

where G^{x,y) is the propagator of the free scalar held, and we hnd immediately that G'{x,y) 
satishes the equation 

dlG'{x,y) =-Jn{x)G°{x,y). (10) 

Unlike G°(x, y), in general G'{x, y) will not depend only on (x—?/), because the source breaks 
the translation invariance of the system. In order to go to momentum space it will therefore 
be necessary to consider the Fourier transform of G'{x,y) with respect to both arguments. We 
dehne G'{p,k) and Jq{p) as follows: 

G'{x,y) = Jd^pJ d^ke^P^e^^y&ip.k) 
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and 

Mx) = Jd^pe^P^Mp), G%x,y) 
The right hand side of m can be rewritten as 


= / d'^k 


^-ik{x-y) 

P ' 


Jn{x)G^{x,y) = d'^p d'^k e'-P^ Jn{p)- 


^-ik(x-y) 


d^k I 


and we obtain the following algebraic equation for the double Fourier transform of the hrst 
order correction to the propagator: 

Transforming back, in conclusion we hnd eq. (Q) of the main text, namely: 

G'{x,y) = I d^P I d^k 

Therefore, if we know the Fourier transform of the characteristic function Jq of the spacetime 
region where the constraint is imposed, we can in principle compute the leading order correction 
to the held propagator and thus to hF[J]. 

It is known 0 that the vacuum-to-vacuum amplitude W[J] = (0'''|0“)j of a held system in 
the presence of an external source J is related to the logarithm of the systems’ ground state 
energy: 

Eo[J] = -T-MnhF[J], 


where the functional integral is supposed to be suitably normalized and the source vanishes 
outside the temporal interval [—T/2,+T/2], with T eventually approaching inhnity. (We use 
units in which % = c = 1.) 

An interesting application of ([TT|) occurs in the case when the held 0(a;) also interacts with 
N static pointlike sources placed at xi,X 2...X7v. Namely, let us add a further, linear coupling 
term Sq to the action of the system: 

N 

Sq = d'^x Q{x)(l){x), with Q{x) = — x,). 

i=i 

The ground state energy of the system corresponds, up to a constant, to the static potential 
energy of the interaction of the sources through the held 0. As before, it is obtained from the 
functional average of the interaction term, computed keeping the constraint into account: 


Eo[J,Q] = t/(xi, ...,X 7 v) =-T Mn(exp {-S'qDj. (12) 
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Expanding ( 0 ) one finds that to leading order in the qjs, f7(xi, is given by a snni of 

propagators integrated on time: 

- ^ f f 
f/(xi, ...,XAr) = -T ^ W / dtj / 

j,i=^ 

where tj, ti e [—T /2, +T/2]. Since the regions Qi are inhnitely elongated in the temporal direc¬ 
tion the fnnction Jq(p -|- k) gets factorized as 

Jnip + k) = (27r)^5(po + ko)jn{p + k). (13) 


Clearly the potential is distnrbed by the presence of the “barriers” jn(x). To hrst order in 
7 we can write 

f/(xi, ...,XAr) = 17°(xi, ...,XAr) +'yU\xi,...,XN) 

and taking into acconnt eq.s (HD, (HD we hnd 

N 


U'{xi,...,xn) = -T ^ W / dtj / dtiG'{xj,xi) = 


il=i 

N 


= -(2'KfT ^ ^ qjQi J dtj J dti J d^p J dk 

il=i 


k2p2 

g*Po(7-ti)+*PXj+*kxiJ^j'p _|_ 


(pg + k2)(p2 + p2^ 


Changing to variables t = tj—ti and s = tj+ti and integrating we hnally obtain the contribntion 
of the pertnrbation to the static potential energy (eq. (|^) of the main text): 


U'{xi,...,xn) = -(27r)® [ dp f dk if^(P + k) 

j,i=i '' 


k2p2 
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e 

P 

Ki 

Pi 

K2 

P2 

K 3 

P 3 

K 4 

P 4 

0.1 

0.3 

12 

10 

80 

10 

12 

10 

80 

10 

0.1 

0.1 

10 

10 

100 

25 

20 

20 

120 

40 

0.1 

0.032 

10 

10 

600 

20 

20 

60 

600 

60 

0.1 

0.01 

20 

20 

1500 

30 

20 

200 

1500 

200 

0.032 

0.032 

10 

10 

2000 

15 

20 

70 

1800 

80 

0.032 

0.01 

12 

12 

6000 

15 

15 

180 

6000 

250 

0.01 

0.032 

10 

10 

7500 

15 

10 

70 

7500 

90 

0.032 

0.0032 

15 

15 

18000 

15 

15 

650 

18000 

700 


Table 1; Integration ranges in the four domains, for some values of e, p. 
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h 

m 

q 

e = 0.1 

(Fig. 1) 

0.14 

0.29 

0.3 

with 

n> 2 

0.12 

0.24 

0.1 

e = 0.05 

(Fig. 3) 

0.17 

0.32 

0.4 

with 

n> 2 

0.16 

0.26 

0.1 


Table 2: Results of the best fit F{n) = exp(—mn + q) + b. 
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FIGURE CAPTIONS. 


Fig. 1 - Dependence of F{e,p,n) on n, in the range n =1...15, for e = 0.1 and p = 0.3. Errors 
are ~ 0.01. 

Fig. 2 - Comparison of the values of F{e,p,n) for e = 0.1 and p = 0.3 (white circles), p = 0.1 
(black triangles) and p = 0.032 (black circles). 

Fig. 3 - Same as in Fig. 1, for e = 0.05. Errors are larger, as shown. 
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